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Abstract 

The weak kaon production off the nucleon induced by neutrinos is studied at the low and in- 
termediate energies of interest for some ongoing and future neutrino oscillation experiments. This 
process is also potentially important for the analysis of proton decay experiments. We develop a 
microscopical model based on the SU(3) chiral Lagrangians. The basic parameters of the model 
are 7,^, the pion decay constant, Cabibbo's angle, the proton and neutron magnetic moments and 
the axial vector coupling constants for the baryons octet, D and F, that are obtained from the 
analysis of the semileptonic decays of neutron and hyperons. The studied mechanisms are the main 
source of kaon production for neutrino energies up to 1.2 to 1.5 GeV for the various channels and 
the cross sections are large enough to be amenable to be measured by experiments such as Minerva 
and T2K. 

PACS numbers: 25.30.Pt,13.15.+g,12.15.-y,12.39.Fe 
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I. INTRODUCTION 



With the recent developments of the atmospheric and accelerator based neutrino experi- 
ments it is now well known that neutrinos oscillate and have finite masses. Now, the main 
goal is to precisely determine the different parameters of the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) matrix, absolute neutrino masses, CP violating phase 6, etc. The neutrino 
energy region of a few GeV is quite sensitive to the neutrino oscillation parameters. There- 
fore, most of the present experiments like MiniBooNE, K2K, T2K, Noz/A, etc. have taken 
data or have been planned in this energy region. Neutrino detection proceeds basically 
through various channels of interaction with hadronic targets like quasielastic scattering, 
meson production, resonance excitations, etc... Therefore, a reliable estimate of these cross 
sections has become important. There is a considerable ongoing theoretical and experi- 
mental effort addressing this question (see e.g. the proceedings of the NUINT conference 
series l|) with many of the studies concentrated at low energies where quasielastic scattering 
and pion production dominate or in the deep inelastic scattering. However, in the discussed 
energy region, other not so well known processes like kaon and hyperon production may also 
become important. In principle, their cross sections are smaller than for the pionic processes 
because of phase space and the Cabibbo suppression for AS' = 1 reactions. Nonetheless, 
in the coming years of precision neutrino physics, their knowledge could be relevant for the 
data analysis, apart from their own intrinsic interest related to the role played by the strange 
quarks in hadronic physics. 

The currently available data is restricted to a few events measured in bubble chamber 
experiments j2-4|. However, this is expected to change soon. In particular, MINERz/A, a 
dedicated experiment to measure neutrino nucleus cross section using several nuclear targets 
like Carbon, Iron and Lead in the neutrino energy region of 1-20 GeV has recently started 
taking data. It is also planned to study specifically the strange particle production and it is 
expected that thousands of events would be accumulated where a kaon is produced in the 
final state 

On the theoretical side there are very few calculations which deal with strange particle 
production at low neutrino energies: single hyperon production 6|, l7|, the study of several 
kaon and hyperon production channels of Dewan [8] and the work of Shrock [9] who has 
analysed the AS = processes. At higher energies, Amer has studied the strange particle 
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production assuming the dominance of s-channel resonant mechanisms [10|. Also, in part 
as a consequence of the scarcity of theoretical work, the MonteCarlo generators used in 
the analysis of the experiments apply models that are not too well suited to describe the 
strangeness production at low energies. For instance, NEUT, used by Super-Kamiokande, 
K2K, SciBooNE and T2K, only considers associated kaon production implemented by a 



model based on the excitation and later decay of resonances ll|] . A similar model is used by 



other event generators like NEUGEN [12|], NUANCE [l3| (see also discussion in Ref. [ij) 
and GENIE jlSs]. As it will be emphasized below, this approach is not appropriate for low 
energies strangeness production. 

In neutrino induced reactions, the first inelastic reaction creating strange quarks is the sin- 
gle kaon production (without accompanying hyperons)^. This charged current (CC) AS = 1 
process is particularly appealing for several reasons. One of them is the important back- 
ground that could produce, due to atmospheric neutrino interactions, in the analysis of one 
le main decay channels the proton has in many SUSY GUT models {p ^ u + K^) 



of t 



16- 



18| . A second reason is its simplicity from a theoretical point of view. At low energies, 
it is possible to obtain model independent predictions using Chiral Perturbation Theory 
(xPT) and due to the absence of S* = 1 baryonic resonances, the range of validity of the 
calculation could be extended to higher energies than for other channels. Furthermore, the 
kaon associated production (with accompanying hyperons) has a higher energy threshold 
(1.10 vs. 0.79 GeV). This implies that even when the associated production is not Cabibbo 
suppressed, for a wide energy region (such as the ANL, the MiniBooNE or the T2K neutrino 
spectrum) single kaon production could still be dominant 

The paper is organized as follows. In Sec. [TTlwe present the formalism for CC single kaon 
production in neutrino nucleon scattering based on the Lagrangians of SU(3) xPT. We also 
discuss the differences with previous calculations. Results, discussions and our concluding 
remarks are presented in Sec. Illli 



^ For antineutrinos the lowest threshold for |A5'| = 1 reactions is much lower and corresponds to hyperon 
production. 
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II. FORMALISM 



The basic reaction for the neutrino induced charged current kaon production is 

vi{k) + N{p) ^ l{k') + N\p') + K{pu): (1) 

where I = e, fi and A^&A^'=n,p. The expression for the differential cross section in the 
laboratory (lab) frame for the above process is given by, 

d'cT = , "^^'j ^k + p-k'-p'- Pk)mM\', 

where k and k' are the 3-momenta of the incoming and outgoing leptons in the lab frame 
with energy E and E' respectively. The kaon lab momentum is pk having energy Ek, 
M is the nucleon mass, SS|A^p is the square of the transition amplitude matrix element 
averaged(summed) over the spins of the initial(final) state. At low energies, this amplitude 
can be written in the usual form as 

M = ^#)J^(^) = -^A'^T^^J''^''\ (2) 
V2 2V2 M^2V2 ^ ^ 

j1 



where j^^-* and J'^(^) are the leptonic and hadronic currents respectively, Gp = a/2 
1.16639(1) X lO^^GeV^^ is the Fermi constant and g is the gauge coupling. The leptonic 
current can be readily obtained from the standard model Lagrangian coupling the W bosons 
to the leptons 
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2^2 l^'^'''^'^^ - ^'^^ + W^h'il - = [j^L)W;^ + h.c.\ . (3) 

We consider four different channels that contribute to the hadronic current. They are 
depicted in Fig. [T] There is a contact term (CT), a kaon pole (KP) term, a u-channel process 
with a S or A hyperon in the intermediate state and finally a meson (vr, 77) exchange term. 
For the specific reactions under consideration, there are not s-channel contributions given 
the absence of 5* = 1 baryonic resonances. The current of the KP term is proportional to q^. 
This implies, after contraction with the leptonic tensor, that the amplitude is proportional 
to the lepton mass and therefore very small. 

The contribution of the different terms can be obtained in a systematic manner using 



xPT. This allows to identify some terms that were missing in the approach of Ref. 



a 



which 
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FIG. 1: Feynman diagrams for the process vN — t- IN' K. First row from left to right: contact term 
(labeled CT in the text), Kaon pole term (KP); second row: u-channel diagram (CrS, CrA) and 
Pion(Eta) in flight (vrP, (r/P) 



only included the u-channel diagrams in the calculation. The lowest-order SU(3) chiral 
Lagrangian describing the pseudoscalar mesons in the presence of an external current is 

f2 f2 



(4) 



where the parameter /^r = 92.4MeV is the pion decay constant, U is the SU(3) representation 
of the meson fields 



U (x) = exp i 



(j){x) 



\ V2K- V2K^ -f^v ) 



(5) 



and D^U is its covariant derivative 



D^U = d^U~tr^U + tUl^. 



(6) 



Here, and correspond to left and right handed currents, that for the CC case are given 
by 



1, = -^{W+T+ + W-T_), 



(7) 
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with the W boson fields and 



^ Vud Vus ^ 


\0 J 



( ^ 

Vud 

V y 



Here, Yij are the elements of the Cabibbo-Kobayashi-Maskawa matrix. The second term of 
the Lagrangian of Eq. that incorporates the explicit breaking of chiral symmetry coming 
from the quark masses |19], is not relevant for our study. 

The lowest-order chiral Lagrangian for the baryon octet in the presence of an external 
current can be written in terms of the SU(3) matrix 

\ 



B 



.J_yo I J_A 



V 

n 



(8) 



as 



(1) 

MB 



Tr 



B {tip - M) b] - ^Tr (Bnd^,, B}) - ^Tr (3^16^, B]) , (9) 



where M denotes the mass of the baryon octet, and the parameters D = 0.804 and F = 0.463 
can be determined from the baryon semileptonic decays [2^ . The covariant derivative of B 
is given by 

D^B = d,B + [r^,B], (10) 

with 



1 



" 2 

where we have introduced = U. Finally, 



(11) 



(12) 



The next order meson baryon Lagrangian contains many new terms (see for instance 



Ref. 



2l|). Their importance for kaon production will be small at low energies and there 



are some uncertainties in the coupling constants. Nonetheless, for consistency with previous 
calculations, we will include the contribution to the weak magnetism coming from the pieces 

C'Sb = 4Tr , a^'-B]) + d.Ti , a^'^B}) + . . . , 

where the tensor can be reduced for our study to 

fjiu d^li, di/lfj^ '^[^/ii ^u\- 
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(13) 



(14) 



TABLE I: Values of the parameters appearing in the hadronic currents. 



Process 


Act Bct ^CrS ^CrA Akp At,p Arjp 


un —7- IKn 


1 D-F -(D-F) 111 


up — )• IKp 


2 -F -(D-F)/2 (D+3F) 2-11 


un — )• IKp 


1 -D-F (D-F)/2 (D+3F) 1 -2 



In this case, the coupHng constants are fully determined by the proton and neutron anoma- 
lous magnetic moments. This same approximation has also been used in calculations of 



single pion production induced by neutrinos 22|]. Now, writing the amplitude for the cou- 
pling of the W boson to the hadrons for each of the terms in the form ^^{Jh^^ + h.c), 
for consistency with Eq. HI we get the following contributions to the hadronic current 

r\c^ = -iActVus^N{p'){Y + ri''BcT)N{p), 



{p - Pkf - Mi 

= iA,pVUD - 3i^)^ (^_^^^_^2 ^(p07^(g^ - 2pfc^)iV(p), (15) 

where, q = k — k' is the four momentum transfer, Vus = sinO = 0.22 where 6 is the Cabibbo 
angle, N{- ), N{- ) denote the nucleon spinors, fip = 1.7928 and yU„ = —1.9130 are the proton 
and neutron anomalous magnetic moments. The value of the various parameters of the 
formulas are shown in Table [H One can notice the induced pseudoscalar form factor in the 
j^\cri:crA curreuts, which takes into account the coupling of the W boson to the baryon 
through a kaon. However, as for the KP term, its contribution is suppressed by a factor 
proportional to the final lepton mass and is negligible. Now, we discuss in some detail the 



terms that appear in the couphng of the weak currents to the octet baryons in the u-channel 
diagrams. With very general symmetry arguments, this couphng can be described in terms 
of three vector and three axial form factors. Following the notation of Ref. Q| we have 



0(t 



B 
92 



a^'qu + 



A 

Mb 

93 



q' 



(16) 
(17) 



Mb " Mb 

where Mb is the baryon mass. At the order considered, the chiral Lagrangian provides finite 
values for /i, the weak magnetism form factor /2, gi and a pole contribution to (73. The scalar 
/a and a non-pole part of the pseudoscalar (73 form factors would only appear at higher orders 
of the chiral expansion. Furthermore, their contribution to the amplitude is suppressed by a 
mi (lepton mass) factor and they are usually neglected. The value of (72 vanishes in the limit 
of exact SU(3) symmetry and there is very little experimental information about it. In fact. 



it is also neglected in most analyses of hyperon phenomenology 



23| . The values of /i and 



gi obtained from the lowest order chiral Lagrangians describe well the hyperon semileptonic 



decays 



20 



23 



24| 



Eventually, if the cross sections for the discussed processes were measured with some 
precision, one could use them to explore these form factors at several values. The current 
experimental information, based on the semileptonic decays, covers only a very reduced 
range for this magnitude. 

Finally, we consider the dependence of the weak current couplings provided by the 
chiral Lagrangians discussed earlier. We should remark that, even at relatively low energies 
and low momenta of the hadrons involved in our study, reaches moderate values. The 
dependences of the needed form factors (e.g. Kvr, YN) are poorly known if at all. Several 
prescriptions have been used in the literature. For instance, for quasielastic scattering and 



single pion production, the vector 
electromagnetic ones (see e.g. 
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orm 
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actors are usually related to the well known nucleon 



26( 1 and references therein). This procedure is well 



suited for these two cases because of isospin symmetry. However, in the SU(3) sector we 
expect to have some symmetry breaking effects. Similarly, for the axial form factors, a 
q^ dependence obtained from the nucleon-nucleon transition obtained in neutrino nucleon 
quasielastic scattering is normally used. However, the axial mass is not well established and 
it runs from values around 1 GeV 
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28| to 1.2 GeV recently obtained by the K2K 



and 



MiniBooNE [30| collaborations. Again here, we expect a different behavior for the hyperon- 




(GeV) 



FIG. 2: Contribution of the different terms to the total cross section for the z/^p — ?■ jiK^p reaction. 



nucleoli vertices. One of the possible choices (e.g. [lO|) is to use a dipole form with the mass 
of the vector(axial) meson that could couple the baryon to the current. In this work, in view 
of the present uncertainties, we adopt a global dipole form factor -F(g^) = 1/(1 — /MpY, 
with a mass Mp — 1 GeV that multiplies the hadronic currents. Its effect, that should be 
small at low neutrino energies will give an idea of the uncertainties of the calculation and 
will be explored in the next section. 



III. RESULTS AND DISCUSSION 

We consider the following reactions: 

ui+p ^ r + K+ +p (/ = e, /i) (18) 
ui + n r + K° +p 
vi + n -> r + + n 

The total scattering cross section a has been obtained by using Eq. ([2]) after integrating 
over the kinematical variables. In Figs. (jSjH]), we present the results of the contributions of 
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FIG. 3: Contribution of the different terms to tlie total cross section for the i/^n — )• jiK^p reaction. 



the different diagrams to the total cross sections. The kaon pole contributions are negligible 
at the studied energies and are not shown in the figures although they are included in the 
full model curves. We observe the relevance of the contact term, not included in previous 
calculations. Starting from the z/^ + p — i- ^~ + + p channel, we find that the contact 
term is in fact dominant, followed by the u-channel diagram with a A intermediate state and 
the TT exchange term. As observed by Dewan Q| the u-channel S contribution is much less 
important, basically because of the larger coupling {NKA ^ NKH) of the strong vertex. 
The curve labeled as Full Model has been calculated with a dipole form factor with a mass of 
1 GeV. The band corresponds to changing up and down this mass by a 10 percent. A similar 
effect is found in the other channels and we will only show the results for the central value 
of 1 GeV. We have also checked that the cross section obtained without the contact term 
and after correcting for the different values of the Cabibbo angle and the Yukawa strong 
coupling agrees well with the result of Fig. 7 of Ref. js] at its lowest energy. Higher energies 
are well beyond the scope of our model. 

The process z/^ + — )■ /i~ + + p has a cross section of a similar size and the contact 
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FIG. 4: Contribution of the different terms to the total cross section for the v^n — )• fiK^n reaction. 

term is also the largest one, followed by the vr exchange diagram and the u-channel (A) term. 
The rate of growth of the latter is somehow larger and could become more important at 
higher energies. As for the previous channel, we observe a destructive interference between 
the different terms and the cross section obtained with the full model is smaller than that 
produced by the contact term alone. 

Finally, the reaction z// + — > Z~ + + n has a smaller cross section. The pion excharige 
term is substantially bigger than the u-channel mechanisms, as already noted in Ref. [8|. 
The contact term is also dominant for this channel and the total cross section calculated 
only with this term practically coincides with the full result. Therefore, we have found that 
the contact terms, required by symmetry, play a major role in the description of the kaon 
production induced by neutrinos at low energies. 

Above the energy threshold for the production of kaons accompanied by hyperons, this 
latter kind of processes could have larger cross sections due to the larger coupling for AS" = 0, 
(Vud vs Vus)- To explore this question and the range of energies where the processes we have 
studied are relevant we compare our results in Fig. |5l with the values for the associated 

n 

production obtained by means of the GENIE Monte Carlo program p^]. We observe that, 
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FIG. 5: Cross sections as a function of the neutrino energy for single kaon production vs. associated 
production obtained with Genie jisl. 



due to the difference between the energy thresholds, single kaon production for the z/^ + p — )■ 
l~ + + p is clearly dominant for neutrinos of energies below 1.5 GeV. For the other 
two channels associated production becomes comparable at lower energies. Still, single 
production off neutrons is larger than the associated production up to 1.3 GeV and even the 
much smaller K'^ production off neutrons is larger than the associated production up to 1.1 
GeV. The consideration of these AS* = 1 channels is therefore important for the description 
of strangeness production for all low energy neutrino spectra and should be incorporated in 
the experimental analysis. 

e [TT] we show the total cross section results for the three channels averaged over 
311 ]. the MiniBooNE [321] and the off-axis (2.5 degrees) T2K [3ii] muon neutrino 



In Tab 
the ANL 



fluxes, all of them peaking at around 0.6 GeV. After normalization of the neutrino flux 
we have 



dE4>{E)(r{E), 



(19) 



where ii^th is the threshold energy for each process and -Ehigh is the maximum neutrino 
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TABLE II: Cross sections averaged over the neutrino flux at different laboratories in units of 10 
cm^. Theoretical uncertainties correspond to a 10% variation of the form factor mass. 



Process 


ANL MiniBooNE T2K 


i/fj^n — )■ K^n 


0.06(1) 0.07(1) 0.09(1) 




0.28(5) 0.32(5) 0.43(8) 


Vfj^n —7- fi~ K^p 


0.17(3) 0.20(3) 0.25(5) 



TABLE III: Number of events calculated for single kaon production in water corresponding to the 
SuperK analysis for atmospheric neutrinos. 



Process 


Events e Events fj. 


f/n — )• l~nK^ 
vin — )■ l^pK^ 
vip —7- l^pK^ 


0.16 0.27 
0.45 0.73 
0.95 1.55 


Total 


1.56 2.55 



energy. As discussed previously, in these three cases, the neutrino energies are low enough 
for single kaon production to be relevant as compared to associated kaon production. Also 
the invariant mass of the hadronic system and the transferred momentum only reach the 
relatively small values where our model is more reliable. 

We can get an idea of the magnitude of these channels by comparing their cross section 
to some recent results. For instance, the cross section for neutral current vr" production per 
nucleon has been measured by the MiniBooNE collaboration 3J] obtaining a = (4.76±0.05± 



0.76) X 10"''° cm^ with a data set of some twenty thousand valid events. The cross sections 
predicted by our model with the same neutrino flux are around two orders of magnitude 
smaller, what means that a few hundreds of kaons should have been produced. 

The atmospheric spectrum [35j also peaks at very low energies and our model should be 
very well suited to analyse the kaon production. In Table llllt we show the number of kaon 
events that we obtain for the 22.5 kTons of a water target and a period of 1489 days as in 



the SuperK analysis 
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36| of proton decay. As in the quoted paper, we include cuts in 



the electron momentum (pe > 100 MeV) and muon momentum (pe > 200 MeV). We find 
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that single kaon production is a very small source o: 



37 



background. In the SuperK analysis 
and only included associated kaon 



the kaon production was modeled following Ref. 
production. Although some of the cuts applied in their analysis, such as looking for an 
accompanying hyperon, are useless for our case, we find that this source of background is 
negligible, given the smallness of our results and the totally different energy distribution of 
kaons and final leptons in the production and decay reactions. 

Finally, we study the values of involved in the reaction for the typical neutrino energies 
we have considered. If high values of this magnitude are relevant, the results would be 
sensitive to higher orders of the chiral Lagrangians and/or a more precise description of the 
form factors. We show the distribution in Fig. (jS]) for the three studied channels at a 
neutrino energy E,^ = 1 GeV. The reactions are always forward peaked (for the final lepton), 
even in the absence of any form factor {F{q^) = 1), favouring relatively small values of the 
momentum transfer. In this figure, we also show the dependence of the cross section on 
the mass of the final lepton that reduces the cross section at low values. The process 
Ue + n ^ + + p shows a slightly different behavior that reflects an important (and 
dependent) interference between the pion exchange and the contact terms. 

Till now we have discussed the kaon production off free nucleons. However, most of the 
experiments are carried out on detectors containing complex nuclei such as iron, oxygen or 



carbon. On the other 
induced by neutrinos 



land, nuclear effects are known to be quite large for pion production 



26l . I39l-l4l|. Fortunately, this question is much simpler for the kaons. 



First, because there is no kaon absorption and the final state interaction is reduced to a 
repulsive potential, small when compared with the typical kaon energies. Second, because 
of the absence of resonant channels in the production processes. We could remember here 
that some of the major nuclear effects for pion production are originated by the modification 
of the A(1232) properties on nuclei. Other nuclear effects, such as Fermi motion and Pauli 
blocking will only produce minor changes on the cross section and can easily be implemented 
in the Monte Carlo codes. 

In summary, we have developed a microscopical model for single kaon production off 
nucleons induced by neutrinos based on the SU(3) chiral Lagrangians. This model should 
be quite reliable at low and intermediate energies given the absence of = 1 baryonic 
resonances in the s-channel. The parameters of the model are well known: /^r, the pion 
decay constant, Cabibbo's angle, the proton and neutron magnetic moments and the axial 
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(GeV^) 

FIG. 6: at = IGeV for single kaon production induced by neutrinos. The curves are labeled 
according to the final state of the process. 

vector coupling constants for the baryons octet, D and F. For the latter ones, we have taken 
the values obtained from the analysis of the hyperon semileptonic decays. The importance 
of higher order terms has been estimated using a dipole form factor with a mass around 
IGeV and exploring the dependence of our results on this parameter. 

We obtain cross sections that are around two orders of magnitude smaller than for pion 
production for neutrino spectra such as those of ANL or MiniBooNE. This can be understood 
because of the Cabibbo suppression and of the smaller phase space. Nonetheless, the cross 
sections are large enough to be measured, for instance, with the expected Minerva and T2K 
fluxes and could have been well measured at MiniBooNE. We have also found, that due to 
the higher threshold of the associated kaon production, the reactions we have studied are 
the dominant source of kaons for a wide range of energies, and thus their study is important 
for some low energy experiments and for the atmospheric neutrino flux. 
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